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ABSTRACT 



Rapid accretion of cold intergalactic gas plays a crucial role in getting gas into galaxies. It has been 
suggested that this gas accretion proceeds along narrow streams that might also directly drive the turbu- 
lence in galactic gas, dynamical disturbances, and bulge formation. In cosmological simulations, however, 
it is impossible to isolate and hence disentangle the effect of cold stream accretion from internal instabili- 
ties and mergers. Moreover, in most current cosmological simulations, the phase structure and turbulence 
in the interstellar medium (ISM) arising from stellar feedback are treated in an approximate (sub-grid) 
manner, so that the feedback cannot generate turbulence in the ISM. In this paper we therefore test the 
effects of cold streams in extremely high-resolution simulations of otherwise isolated galaxy disks us- 
ing detailed models for star formation and stellar feedback; we then include or exclude mock cold flows 
falling onto the galaxies, with mass accretion rates, velocities and flow geometry set to maximize their 
effect on the gaseous disk. We find: (1) Turbulent velocity dispersions in gas disks are identical with or 
without the presence of the cold flow; the energy injected by the flow is efficiently dissipated where it 
meets the disk. (2) In runs without stellar feedback, the presence of a cold flow has essentially no effect 
on runaway fragmentation (local collapse), resulting in star formation rates (SFRs) that are an order-of- 
magnitude too large. (3) Model disks in runs with both realistic feedback and cold flows have higher SFRs, 
but only insofar as they have more gas. (4) Because the flows are extended, relative to the size of the disk, 
they do not trigger strong resonant responses and so induce weak gross morphological perturbation (bulge 
formation via instabilities/fragmentation is not accelerated). (5) However, flows can thicken the disk by 
direct contribution of out-of-plane or mis-aligned star-forming streams/filaments. We conclude that while 
inflows are critical over cosmological timescales to determine the supply and angular momentum of gas 
disks, they have weak instantaneous dynamical effects on galaxies. 
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1 INTRODUCTION 



Since some of the earliest models of galaxy formation in a cos- 
mological setting, it has been known that in low-mass halos (< 
10 11-12 Mq), the gas cooling time is shorter than the dynamical 
time even if gas is shock heated to virial temperature before accre- 
tion. In this case shock heating during virialization will not provide 
sufficient pressure support for the infalling gas and gas could reach 
galaxies in near free-fall (|Silk||1977[ |Rees & Ostriker|1977| [BnTl 
|ney| 1977] | White & Rees|1978[ for a more detailed spherical shock 
stability analysis with a modern ACDM cosmology, see Birnboim 
& Dekel 2003 ). In recent years, numerical hydrodynamic simula- 
tions have considerably improved our understanding of this regime 
which is part of a more general process termed "cold accretion" or 
"cold flows." Keres et al. (2005 ) argued that in lower mass sys- 
tems, gas can reach galaxies in free-falling filaments which are 
cold (~ 10 4 — 10 5 K) by virtue of their large densities preventing 
an accretion shock; at high redshifts (z > 2), such streams may 
even penetrate more massive halos (where the spherically- averaged 
^cooi ^> fdyn) because of the highly non-spherical nature of accretion. 
In this early work the gas accreted by galaxies in massive halos was 
shown to largely come from the "hot mode accretion", i.e. most 
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of the infalling filamentary gas heats before reaching galaxies. A 
large body of subsequent theoretical work has developed on this 
topic (see e.g. |Dekel & Birnboim|[2006l |Dekel et al.||2009b[ [ Op^] 
Ipenheimer et al.|2010[ |Ocvirk et al.|2008[|Keres et ad.|2009| and 
references therein). At the same time, tantalizing observational ev- 
idence for the existence of infalling cool gas is emerging (Ri baudo| 
|et al.| 2011 ; Kacprzak et al. 2012 ), although systematic detection of 
large gas infall rates in high-redshift halos is still lacking ( |Sancisi| 
|et al.|2008l|Steidel et al.|20T0] ). 

The exact nature of the accretion process in massive (~ 
10 12 Mq) halos is still a subject of a debate, owing to the sensitiv- 
ity of predictions to details of numerical methods. All simulations 
agree that most of the gas in massive halos heats to high temper- 
atures. The debate is mostly about the nature of the infalling gas 
in the vicinity of galaxies. At z ~ 2 cold streams of gas can still 
bring material deeply into the halos, however most of this gas is 
heated before it is accreted to the halo center ( [Keres et"aT] 2005, 
|Ocvirk et al.|2008| >. Subsequent work ( [Keres et al.|2009| > using dif- 
ferent hydrodynamic codes found that a only relatively small frac- 
tion (typically ~ 10 — 20% of the total halo infall at z ~ 2, see 
|Faucher-Giguere et al.||201l) survives all the way to the central 
galaxy. Cooling of the hot halo in these simulations was not ef- 
ficient, so this small fraction of surviving cold gas was found to 
dominate the gas supply to the halo center. Using another tech- 
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nique, |Dekel et al.| ( |2009b| ) suggested that at z ~ 2 — 3 rapid infall 
continues unimpeded via narrow filaments all the way to the halo 
center even for massive halos. 

More recent simulations using the moving mesh code AREPO 
(which resolves some numerical inaccuracies in previous work; 
see|Springel 2010 ) find that (absent galactic outflows) the infalling 
gas is indeed filamentary within the virial radius but gets strongly 
heated before accretion onto the galaxy (Nelson et al. 2013). In 
these runs, at z ~ 2 a large fraction of gas in a ~ 10 Mq halo flows 
as cool gas down to ~ 0.3/^; the infalling gas then shock-heats, 
but remains in a coherent structure before it cools again rapidly 
close to the galaxy (see Keres et al.1|2012| |Nelson et al]|2013| ). 
However, the gas heating removes some kinetic energy of the in- 
fall and forms structures whose coherence lengths exceed the size 
of the galactic disk (~ lOkpc). Further complexity comes from the 
halo mass dependence; in lower-mass halos the width of filaments 
is typically large compared to the galaxy (even when filaments are 
cold). Thus, the simplest picture of narrow streams of cold gas im- 
pacting the galactic disk may need to be modified (see discussions 
in |Katz et al. [20031 |Keres et al.|2009| [Nelson et al.|2013) . 

In any case, the consequences of narrow filamentary inflows 

- assuming they exist - are controversial. For example, it is well- 
known that in a gaseous disk where cooling is efficient, turbulent 
velocity dispersions should dissipate in a crossing time, so some- 
thing must continuously maintain the turbulence in galaxies. Al- 
though this is often attributed to a combination of stellar feedback 
and local gravitational instabilities, some authors have argued that 
the observed turbulent velocity dispersions - especially in high- 
redshift galaxies - may be directly driven by the accretion energy 
(e.g. Krumholz & Burkert 2010 ; Cacciato et al. 201 1 ). By extension 
(if the turbulence can be both driven and efficiently isotropized) 
this could argue for turbulence driving gas disk thicknesses/scale- 
heights. In some versions of this argument, maintaining the large- 
scale turbulence would in turn be sufficient to explain the low star 
formation efficiencies of galaxies (i.e. the Kennicutt 1998 relation). 
Others have argued that the flows directly trigger or enhance sec- 
ular instabilities in the disks, driving bulge formation (Deke l et al.| 
|2009bl >. 

If true in an instantaneous sense, this would imply a radical 
departure from one of the fundamental assumptions in most ana- 
lytic or semi-analytic galaxy formation models. In such models, ac- 
cretion is the "fuel supply" which determines the gas mass, angular 
momentum content, radius, etc. of the galaxy disk; these properties 
in turn determine things like the star formation rate and local sta- 
bility properties, hence ultimately the disk structure and feedback 
properties. But the inflowing gas, excluding discrete merger events, 
does not directly alter the structural properties of the disk. If the 
inflow powers the turbulence in the disk, for otherwise identical 
properties in that disk, then the instantaneous properties of the "ac- 
cretion flow" will significantly change the gas velocity dispersions 
within the disk, and hence all the disk properties that depend on the 
velocity dispersion (e.g., the star formation rate and stability). 

It is critical to assess this claim in numerical simulations, 
which can simultaneously follow the relevant physics: inflow, gas 
cooling, star formation, stellar feedback and the formation of ISM 
phase structure. Such simulations can follow the complex, non- 
linear dynamics of turbulence, torques, and energy and angular mo- 
mentum exchange in gas+stellar+dark matter disk galaxies. How- 
ever, in traditional cosmological simulations - which can self- 
consistently follow the formation of inflows and their properties 

- one cannot simply "remove" the accretion flows, nor vary their 
properties systematically while keeping the galaxy properties fixed. 



In other words, cosmological simulations are not controlled exper- 
iments as far as accretion flows are concerned. 

In addition, achieving the resolution needed to resolve the in- 
ternal structural properties of the gas disk (resolving the vertical 
turbulent cascade, for example) is very challenging in a cosmologi- 
cal simulation, with simulation volumes measured in tens of mega- 
parsecs or larger. Finally, it is well-known that the gas properties 
are strongly affected by feedback from massive stars, so including 
realistic models for these feedback properties is critical. 

In this paper, we therefore use a suite of high resolution hy- 
drodynamic simulations of isolated disk galaxies (with properties 
chosen to match observations) and realistic stellar feedback mod- 
els, including various simplified "cold flow" models of infalling 
gas, to investigate these questions. In order to study the influence 
of infalling gas on the gas dynamics in galaxies, we will consider 
the most optimistic assumptions (ignoring some of the subtleties 
above). We will assume that narrow filaments of cool gas do man- 
age to survive and plunge into the disk of the galaxy. Furthermore, 
we assume that infall rates from cold mode accretion onto galaxies 
represent the full infall rate onto halos (clearly an upper limit). Our 
purpose is to see if such filamentary accretion can drive high veloc- 
ity dispersion and scale-heights of galaxies, and to assess the rela- 
tive importance of this process when compared to star formation- 
driven feedback. 

2 THE SIMULATIONS 
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The simulation method we use is described in detail in Hopkins 
< 201 lbj (§ 2 & Tables 1-3; hereafter Paper I) and Hopkins 
(§ 2; hereafter Paper II). We briefly summarize the 
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most important properties here. The simulations were performed 
with a heavily modified version of the parallel TreeSPH code 
GADGET-3 ([Spr ingel 2005 ), in a fully conservative formulation 
(Springel & Hernquist 2002) which is also density-independent 
in a manner that allows contact discontinuities and improved fluid 
mixing ( Hopkins 2012a ). The artificial viscosity, adaptive timestep- 
ping, and smoothing kernel are updated as described therein^They 
include stars, dark matter, and gas, with cooling, shocks, star for- 
mation, and stellar feedback. 

In our simulations gas follows an atomic cooling curve with 
additional fine-structure cooling to 10 K, allowing for the photo- 
ionizing background computed in Faucher-Gigue re et al.| {"2008) 
and gas self- shielding. Metal-line cooling is followed species-by- 
species for 11 tracked species as in |Wiersma et al. ( 2009a. b). The 
enrichment for each species is followed with the appropriate time 
dependence with independent yields directly attached to mass and 
energy return from individual feedback sources described below. 
Star formation is allowed only in dense, molec ular, self-gravitating 
regions above n > 1000 cm -3 . We follow Krumholz & Gnedin 



1 There has been considerable discussion in the literature regarding subtle 
numerical effects on inflow properties and the treatment of sub-sonic mul- 
tiphase flows (see references in §[TJ. Our method is specifically designed 
to better treat this regime. However, we have also re-run our HiZ model 
with feedback with the alternative SPH formulation in Springel & Hern- 
quist ( 2002 1 and that in Abel ( 201 1 ), as well as varied artificial viscosity and 
resolution, and find our major conclusions are robust. We believe this is be- 
cause we focus on the galactic disk, not flow formation and interaction with 
the hot halo. In the former regime, the shocks and relevant turbulent cascade 
are highly super-sonic and external forcing (from gravity and feedback) is 
important; this is the regime where the various numerical approaches agree 
well (see |Kitsionas et al.|2009l|Price & Federrath|2010||Bauer & Spr ingel 
[20121. 
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HiZ No Flow +Flow (Prograde) +Flow (Polar) 




Figure 1. Morphology of the gas in simulations of a massive z ~ 2 — 4 
starburst disk with M* ~ 200 Mq yr _1 . Images are shown at ~ 3.5 orbital 
times when the disk is in a feedback-regulated steady-state. We show face- 
on (upper) & edge-on (lower) projections, at large (halo) scales & small 
(star-forming disk) scales, as labeled. All models shown include stellar 
feedback. We compare a case with no inflow (left), with our "standard" 
mock cold flow (inflow rate M ~ 200 Mq yr _1 , aligned with disk rotation; 
center), and with a similar flow oriented perpendicular to the disk rotation 
(right). The images are at extremely high resolution/depth. Brightness en- 
codes projected gas density (logarithmically scaled with w 6dex stretch); 
color encodes temperature with blue/white being T < 1000K molecular 
gas, pink ~ 10 4 - 10 5 K warm ionized gas, and yellow > 10 6 K hot gas. 
The qualitative morphology, with massive kpc-scale molecular complexes, 
is similar in each case. Warps and out-of -plane streams make the disk ap- 
pear thicker with flows. Out-flowing winds are seen in all simulations, but 
are slightly enhanced when inflows are present. 

( |2011| > to calculate the molecular fraction fn 2 in dense gas as a 
function of local column density and metallicity, and allow SF only 
from molecular gas. We also restrict star formation to gas which 
is locally self-gravitating, i.e. has a = 5v 2 5r/Gm gas (< 6r) < 1 on 
the smallest available scale (Sr being our force softening or smooth- 
ing length). This forms stars at a rate = /w/fe (i.e. 100% ef- 
ficiency per free-fall time); however the galaxy- wide efficiency is 
generally much lower because of feedback. In Paper I and Paper II 
we show the galaxy structure and SFR are basically independent of 
the small-scale SF law, density threshold (provided it is high), and 
treatment of molecular chemistry, because the star formation (and 
hence galactic structure) is feedback-regulated. 

Stellar feedback is included, from a variety of mechanisms: 
(1) Local Momentum Flux from radiation pressure, super- 
novae, & stellar winds: Gas within a GMC (identified with an on- 
the-fly friends-of-friends algorithm) receives a direct momentum 
flux from the stars in that cluster/clump. The momentum flux is 



P — ^SNe + Pw + ^Prad, where the separate terms represent the di- 
rect momentum flux of SNe ejecta, stellar winds, and radiation 
pressure. The first two are directly tabulated for a single SSP as 
a function of age and metallicity Z and the flux is directed away 
from the stellar center. Because this is interior to clouds, the sys- 
tems are always optically thick, so the latter is approximately 

Prad « (1 + Ti R )Lincident/c, where 1 + TjR = 1 + Sgas KjR aCCOUntS 

for the absorption of the initial UV/optical flux and multiple scat- 
terings of the IR flux if the region is optically thick in the IR (with 
Sgas calculated for each particle). 

(2) Supernova Shock-Heating: Gas shocked by supernovae 
can be heated to high temperatures. We tabulate the SNe Type-I and 
Type-II rates from |Mannucci efaT] ( [2006] ) and STARBURST99, 
respectively, as a function of age and metallicity for all star particles 
and stochastically determine at each timestep if a SNe occurs. If so, 
the appropriate mechanical luminosity is injected as thermal energy 
in the gas within the nearest ~ 32 neighbors of the star particle. 

(3) Gas Recycling and Shock-Heating in Stellar Winds: 
Gas mass is returned to the ISM from stellar evolution, at a rate tab- 
ulated from SNe and stellar mass loss (integrated fraction « 0.3). 
The SNe heating is described above. Similarly, stellar winds are 
assumed to shock locally and inject the appropriate tabulated me- 
chanical luminosity L(t, Z) as a function of age and metallicity into 
the gas within a smoothing length. 

(4) Photo-Ionization & Photo-Electric Heating: We also 
tabulate the rate of production of ionizing photons for each star 
particle; moving radially outwards from the star, we then ionize 
each neutral gas particle (using its density and state to determine 
the necessary photon number) until the photon budget is exhausted. 
Ionized gas is maintained at a minimum ~ 10 4 K until it falls out- 
side an HII region. Photo-electric heating is followed in a similar 
manner using the heating rates from Wolfire et al. (1995). 

(5) Long-Range Radiation Pressure: Photons which escape 
the local GMC (not accounted for in (1)) can be absorbed at larger 
radii. Knowing the intrinsic SED of each star particle, we attenu- 
ate integrating the local gas density and gradients to convergence. 
The resulting "escaped" SED gives a flux that propagates to large 
distances, and can be treated in the same manner as the gravity tree 
to give the local net incident flux on a gas particle. The local ab- 
sorption is then calculated integrating over a frequency-dependent 
opacity that scales with metallicity, and the radiation pressure force 
is imparted. 

Details and numerical tests of these models are discussed in 
Paper II. All energy, mass, and momentum-injection rates are taken 
as-is from the stellar population models in STARBURST99, as- 
suming a |Kroupa|{2002) IMF, without any free parameters. Sub- 
tle variations in the implementation do not make significant differ- 
ences to our conclusions. Most important, we do not "turn off" or 
otherwise alter any of the cooling or hydrodynamics of the gas. 

We implement the model in three distinct initial disk mod- 
els spanning a range of galaxy types. Each has a bulge, stellar and 
gaseous disk, halo, and central BH (although to isolate the role of 
stellar feedback, models for BH growth and feedback are disabled). 
At our standard resolution, each model has ^0.3 — 1 x 10 8 total 
particles, giving particle masses of 500 — 1000 Mq and 1 — 5pc 
smoothing lengths, and are run for a few orbital times each. In 
Paper II, a couple ultra-high resolution runs for convergence tests 
employ ~ 10 9 particles with sub-pc resolution. The disk models 
include: 

(1) SMC: an SMC-like dwarf, with baryonic mass Mbar = 
8.9 x 10 8 M (gas m g = 7.5 x 1O 8 M , bulge M b = 1.3 x 1O 8 M , 
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the remainder in a stellar disk m<i) and halo mass Mhaio = 2 x 
1O 1O M . The gas (stellar) scale length is h 8 = 2.1 kpc (h d = 0.7). 

(2) MW: a MW-like galaxy, with halo M ha i = 1.5 X 10 12 , and 
baryonic (M bar , m b , m d , m g ) = (7.1, 1.5, 4.7, 0.9) x 1O 1O M with 
scale-lengths (/y, = (3.0, 6.0) kpc. 

(3) HiZ: a high-redshift massive starburst disk, typ- 
ical of massive star-forming galaxies at z ~ 2 — 4; 
Mhaio = 1.4 x 10 12 Mq (scaled for z = 2 halos), and baryonic 
(Mbar, mb, rrid, m g ) = (17,7,3,7) x 10 10 Mq with scale-lengths 
(fc,^) = (1.6,3.2)kpc. 

We consider each model with and without a "cold flow." Mo- 
tivated by the cosmological simulations cited in § [T] we initialize 
the flow as a narrow cone of gas with opening angle 0, infalling at 
the free-fall velocity (from infinity). We typically assign the flow a 
small initial turbulent dispersion, « 20% of the mean velocity, and 
a sub-virial initial T « 10 5 K; we have varied these parameters, and 
found little effect on the results. The flow is initialized so the mean 
M through all conical annuli is constant, towards an impact param- 
eter at a radius b in the disk, with a relative inclination 6. We ini- 
tialize the flow with sufficient mass and extent to ensure it is main- 
tained for the duration of the simulation. Because such flows are 
believed to be the ultimate source of the disk gas, our default choice 
for b is R e , g as (the effective radius of the gas disk), with a prograde 
orientation (i.e. with the flow having the same sign of angular mo- 
mentum as the disk), with a narrow (20°) opening angle, and M ap- 
proximately equal to the maximal halo baryon accretion rate (halo 
growth rate, for a mean halo of the given mass and redshift, times 
the universal baryon fraction; = 0.15, 10.0, 200MQyr _1 for the 
SMC, MW, HiZ models, respectively). As explained in §[T] these 
rates are purposely set to be maximally optimistic; typical (median) 
halos will have lower accretion rates, and not all of the infall will 
survive all the way to the galaxy as part of a coherent structure. 

3 RESULTS 

Fig. [T] compares the gas morphologies of runs with and without 
"flows" present. The results are broadly similar within the effective 
radii of the gaseous disks. Inflows clearly lead to more gas dis- 
tributed outside the star-forming disk, both in the form of streams 
and filaments falling in, as well as in additional wind material. Over 
several orbits the inflows modify the outer disk, essentially by re- 
placing it with material on new orbits (determined by the flow pa- 
rameters). 

In Figs. |2|3[ we examine the effects on galaxy properties. For 
each model, we run simulations with and without stellar feedback, 
and with and without inflows. In Fig. [5] we vary the properties 
of the flow. As shown in Paper I & Paper II, without feedback, 
the gas radiates away its support and locally collapses into dense 
clumps, globally turning most of its mass into stars in a single dy- 
namical time. The SFRs correspond to M* ~ M gas /%n, an order- 
of-magnitude or more in excess of the observed Kennicutt relation. 
The absolute SFRs decline at later times only because the gas is ex- 
hausted, but always exceed the rate given by the Kennicutt relation. 
The presence of an accretion flow makes no difference to the star 
formation rate relative to Kennicutt; it simply provides more fuel to 
sustain the runaway SFR longer. With feedback present, the SFRs 
self-regulate at a level in good agreement with the observed Ken- 
nicutt relation (Paper I). When feedback is included, the inflow has 
no effect on the SFR, scaled to Kennicutt; at late times the flows 
maintain the disk gas densities, slowing the decline in star forma- 
tion rate. 

We next examine the internal gas disk properties. The three- 
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Figure 2. Effects of cold "accretion flows" on galaxy properties. Each 
column shows a different galaxy type (§ |2j. For each, we show models 
with/without a flow and with/without stellar feedback. Top: SFR vs. time. 
With feedback, models self -regulate near observed SFRs; late-time SFRs 
with flows are higher because gas is replenished, but both lie on the same 
Kennicutt- Schmidt relation. Without feedback, both experience runaway lo- 
cal collapse/SF. Second: Vertical disk scale height h/R. Feedback maintains 
higher values by powering turbulence and isotropizing in-plane gas mo- 
tions. The flows may help isotropise the gas (especially without feedback) 
and boost h /R by depositing low-binding energy material above/below the 
disk. Third: Gas velocity dispersion vs. radius (time- averaged). All mod- 
els maintain Q ~ 1, giving similar a, independent of feedback and flowsP 
Bottom: Final stellar mass profile (at the end of the simulation). Disk in- 
stabilities drive some bulge growth in each case, sensitive to feedback, but 
without a clear dependence on flows. 



dimensional gas velocity dispersion is plotted (as a function of 
radius) in the third row of Fig. [2] We take the average over all 
snapshots after the first disk orbital period, when the system is in 
equilibrium. In Paper Ilwe showed that disks always self-regulate 
at Q ~ 1, so the dispersion a is similar regardless of the pres- 
ence/absence or details of feedback. Here, we see that the presence 
of accretion flows makes almost no difference to the internal veloc- 
ity dispersions of the disks when compared to the run with feedback 
but lacking an inflow]^] 

The second panel of Fig. [5] shows the disk scale height h/R, 
which is related specifically to the vertical disk dispersion a z . In 
Paper II, we showed that h is different in cases with or without 



2 The dispersion in the HiZ case with no feedback or inflow is very large 
at < 1 kpc. This is because the extremely rapid star formation (in the ab- 
sence of feedback), without continuous inflow, rapidly exhausts the central 
gas, and the dispersion is dominated by hot (volume-filling) gas and the 
occasional dense gas clump inspiraling near the circular velocity. 
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Figure 3. Galaxy properties as Fig. [2] Here, we focus on the HiZ model 
(with the most extreme inflow) but vary the inflow properties. All models 
include "cold flows." We compare our standard cold flow model (prograde 
- i.e. aligned with the disk, with M gas = 200 Mq yr _1 and impact parame- 
ter b = R e ) to models with different relative flow/disk orientations (a model 
where the flow is polar, i.e. perpendicular to the disk rotation, and a model 
where the flow is retrograde, i.e. directly oppose the disk rotation), with 
different flow rates (prograde, 400 Mq yr _1 ), and with different impact pa- 
rameters ("narrow" with b = 0.5 R e and cross-section smaller by 1/2, and 
"high- j/wide" with b = 2R e and cross section larger by 2). There is very lit- 
tle dependence of the disk properties on these choices except, perhaps, just 
at the location where the flow meets the disk. Flows with higher angular 
momentum (impact parameter > R e ) are similar but appear to have weaker 
effects on scale-height. 



feedback. Simulations lacking feedback maintain a via global grav- 
itational instabilities and "one-way" collapse, as, e.g. a collapsing 
clump increases in a as it collapses to a resolution-limited small 
radius. In a disk this gravitationally induced velocity dispersion is 
preferentially in-plane and leads to large velocity anisotropy, and 
hence smaller h/R than seen in observations. In the simulations in- 
cluding inflows, the difference between feedback and no-feedback 
runs is reduced somewhat, i.e., the red dotted line recovers to the 
solid red line with a couple orbits. Visually (Fig. [TJ, the disks ap- 
pear "thicker" because the streams directly contribute some gas in 
out-of-plane orbits, and mis-alignment of disk and stream leads 
over a few orbits to warps; the flows may also help isotropise the 
gas dispersions, boosting h/R over what it would be without in- 
flows. Comparing to the run lacking both feedback and inflows, 
the difference is significant. However, these inflow-driven effects 
do not much change h/R compared to the no-inflow run with feed- 
back. And the effect on h/R becomes weaker if we systematically 
increase the impact parameter (and angular momentum) of the in- 
flow; since the dispersions are similar, this appears to be related to 
the ability of the disk to relax to equilibrium. 

We also show the stellar mass profiles at the conclusion of 
each simulation, after ~ 0.5 Gyr of evolution. In Hop kins et al.| 
( |2011a| >, we discuss how feedback alters the efficiency of secular 
evolution, particularly for the HiZ model, comparing the mass pro- 
files with and without feedback. Without feedback, gas experiences 
runaway local collapse into dense "nuggets" which then have no 
choice but to sink to the galaxy center. This leads to artificially 
large/rapid bulge formation. We see here the large difference in the 
bulge formed in the HiZ model with and without feedback (at small 
radii in the bottom row, the black dotted line lies well above the 
other lines). Also, without feedback, the outer disk evolves outward 
as a result of absorbing the angular momentum from these sinking 
clumps. It has been argued that inflows lead to rapid bulge build up 
Dekel et al. (2009a). We do not see this in our simulations — in fact 
if anything we find that inflow, in the absence of feedback, leads 
to smaller bulges than in simulations lacking both feedback and in- 
flows. If feedback is included, the presence or absence of an inflow 
has little or no effect on bulge size. 

4 DISCUSSION 

We consider idealized simulations of isolated disk galaxies with 
and without large "cold streams", where the cosmological inflows 
are represented as simplified streams of cold gas free-falling onto 
the galaxy. This is a crude approximation to the full cosmological 
case, but unlike in cosmological simulations, our simulations allow 
us to add, remove or modify flows while maintaining otherwise 
identical galaxy properties, i.e., they allow for controlled numerical 
experiments. The smaller simulation box sizes used also allow us to 
reach pc- scale and ~ lOOyr resolution and include explicit, spatial 
and time-resolved models for stellar feedback from SNe (Types I & 
II), stellar winds (O-star and AGB), photoionization, and radiation 
pressure. 

We show that in an instantaneous sense - i.e. for otherwise 
fixed galaxy properties - the presence or absence of an inflow has 
only a small effect on the galaxy. The SFR in either case is reg- 
ulated by feedback: as shown in Paper II, with realistic feedback, 
gas cools and collapses until enough young stars are formed to in- 
ject sufficient momentum and regulate against further collapse. As 
a result, the SFR is set by feedback entirely independent of e.g. the 
small-scale star formation law or the presence of cold stream ac- 
cretion. In other words, the flows have no direct effect on the SFR. 
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They do affect the SFR over long times, not because they affect the 
dynamics of the ISM, but because they affect the amount of gas in 
the disk. Without feedback, local collapse runs away in a dynami- 
cal time, giving M* ~ Mg as /%n (°° 50 times larger than observed). 
This is not altered with flows present. In the absence of feedback, 
cold streams are unable to prevent runaway local collapse or "slow 
down" star formation via turbulence. Indeed, as shown in Hopkins 
{20 12b) , any mechanism that drives turbulence only on large scales, 
but does not break up locally self-gravitating regions, is insufficient 
to stave off runaway local collapse. 

We also find that the velocity dispersions are weakly sensitive 
to the flow properties. This is consistent with the result in Paper II 
that a combination of feedback and gravitational instability always 
drives the velocity dispersion to values such that the Toomre Q ~ 1 
(with any "excess energy" efficiently radiated), such that the veloc- 
ity dispersion does not change upon adding or removing individual 
feedback coupling or turbulent driving mechanisms. It follows that 
the stability properties of the galaxies are not significantly altered 
by the presence of inflows. 

Nor do inflows appear to strongly alter the rate of secular evo- 
lution or angular momentum transfer in the disk. Disk instabilities, 
which are present in all the simulations shown here, are primarily 
sensitive to local properties. The inflow we model is smooth and ex- 
tended, and so does not introduce a strong resonant response (see 
e.g.|D'Onghia et al. 2010 ), unlike the case in a galaxy merger. 

Feedback efficiently powers turbulence and isotropises the 
dispersions, so with feedback on we see little difference in the 
disk scale-height when in inflow is turned on or off. However, 
in runs without feedback, the velocity dispersion is largely in- 
plane in the no-flow case, so adding an inflow "stirs" the disk and 
boosts its thickness significantly. We note, though, that these no- 
feedback, inflow- stirred disks are still thinner than disks with feed- 
back present. Flows do make the disks appear thicker, by contribut- 
ing some out-of-plane star-forming gas directly, and, over time, 
generating warps and misalignments. This is not a consequence of 
the flow dynamically disturbing, stirring, or driving the disk, but 
rather of the disk material (coming from the flow) having a dif- 
ferent angular momentum distribution (with mis-aligned material) 
in our idealized, isolated galaxy models. With sufficient time (if the 
flow is stable), the situation will likely relax into a new disk aligned 
with the flow. 

If we examine the location where inflow meets disk, we see 
a mix of behaviors. Some diffuse material undergoes a series of 
radiative shocks, but the relative velocities in each are sufficiently 
small that this places material near the peak of the cooling curve 
and so the energy is efficiently dissipated. More dense streams and 
clumps fall in ballistically, and can penetrate the disk, forming the 
mis-aligned streams discussed above, which are gradually torqued 
into equilibrium orbits. 

We also note that the feedback-driven outflows or winds dis- 
cussed in detail in Hopkins et al. (2012a), are not strongly altered 
with cold flows present, although overall the outflow mass-loading 
rates (M W i n d/M*) are enhanced by factors ~ 1.5 — 2 (reflected in 
the extended outflows in Fig. [TJ, because the inflow includes ma- 
terial which is marginally bound. This marginally bound inflowing 
material can be accelerated by radiation pressure and hot gas ram 
pressure at large radii, i.e., it is accelerated off of the surface of the 
stream. 

The "cold flow" models here represent extreme cosmological 
situations, in that they have mass inflow rates equal to the entire 
baryonic halo accretion rate and assume 100% of the infall is in 
"cold mode." They are inf ailing at the free-fall velocity, reaching 



the galaxy and there is (by design) no hot halo to resist their in- 
fall or strip the streams. They also have impact parameters very 
close to the disk effective radii. Realistic situations in cosmologi- 
cal simulations may produce weaker effects than those here. This 
is especially the case in lower-redshift or more massive halos like 
the MW, where much of the gas is probably shocked to virial tem- 
peratures with long cooling times. Furthermore, in massive halos, 
even at z ~ 2 gas delivered from minor mergers might be impor- 
tant channel of gas supply (Keres et al. 2009 ). Moreover, the recent 
code cosmological simulations mentioned in the Introduction have 
questioned certain previous conclusions about the dynamics of the 
inf ailing gas, arguing for example that much of the angular momen- 
tum transfer associated with infalling gas and supply of the galac- 
tic gas was likely affecte by numerical inaccura cies (fVo gelsberger 
|et al.|2012||Sijacki et al.|20T2l|Keres et al.|20121|BaueT& Springel 
2012 ). These studies have generally found that the streams are more 
smooth and diffuse, and more efficiently stripped/mixed by diffuse 
halo gas as they approach the disk, than previous calculations had 
estimated. 

However, one way in which accretion could introduce stronger 
perturbations than what we consider here is if it were very 
"clumpy," with a large fraction of the incoming mass in units with 
sizeable masses relative to the disk. Hopkins et al. (2008 ) discuss 
the resonant response to radially infalling clumps/galaxies, and 
|D'Onghia et al.|p010| ) consider non-merging passages; both sug- 
gest that significant effects could occur if the clumps have masses 
> 10% of the disk, effectively acting as minor mergers. And in- 
deed there may be a large number of such clumps and/or mergers 
"carried in" with the flows (see e.g. Keres et al. 2005 2009 ; Dekel 
|et al.|2009b| although also see the numerical caveats above). The re- 
sponse of disks to minor mergers and harassment, however, is well- 
developed and studied in a number of other contexts, and is outside 
the scope of our investigation here. But to the extent that this is 
important, our results suggest it may be separable into a merger- 
like contribution from individual massive clumps falling in, and a 
weaker contribution from the "smooth accretion flow." 

We are of course not arguing that "cold" accretion is unimpor- 
tant for galaxy formation. None of our arguments change whether 
gas enters the halo without undergoing virial shocks or with cooling 
times short compared to dynamical times and so reaches the galaxy 
rapidly. This gas is the supply for star formation on cosmological 
timescales; therefore it (with the effects of subsequent feedback) 
determines the disk gas mass, stellar mass, and angular momen- 
tum content of the galaxy, which in turn determine the SFR, veloc- 
ity dispersion, and secular stability. Both theoretical ( [Keres etah] 
2005 ) and observational (e.g. Daddi et al. 2007 ) arguments indi- 
cate that continuous gas accretion is necessary to to explain long 
term star formation activity of galaxies. For example, on cosmo- 
logical timescales, the "equilibrium SFR" of galaxies will be set in 
part by the rate of inflow: the gas consumption rate is M* +M W i n d, 
where M W i n d oc M* for otherwise similar properties (Paper III). If 
this falls below the infall M aC c, a gas reservoir inevitably builds up, 
increasing the SFR until M acc ~ M* +M W i n d- However, SF (and the 
generation of outflows) does not directly "know" about the inflow. 
Rather, inflows (in competition with SF and outflow) determine the 
gas supply, which determines the SFR locally needed to prevent 
runaway local collapse, which determines the actual SFR, which 
then determines the outflow properties. 

In short, we argue that the role of accretion is to cosmolog- 
ically determine the "initial" global properties of the galaxy, e.g., 
the gas and angular momentum supply rate. For a given set of those 
global properties, the local instantaneous properties of the galaxy 
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are set by a combination of local dynamical and feedback processes 
that operate on much smaller spatial and time scales. 
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